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ABSTRACT 

HR8799 is a A Bootis, 7 Doradus star hosting a planetary system and a debris disk 
with two rings. This makes this system a very interesting target for asteroseismic 
studies. This work is devoted to the determination of the internal metallicity of this 
star, linked with its A Bootis nature (i.e., solar surface abundances of light elements, 
and subsolar surface abundances of heavy elements), taking advantage of its 7 Do- 
radus pulsations. This is the most accurate way to obtain this information, and this 
is the first time such a study is performed for a planetary-systeni-host star. We have 
used the equilibrium code CESAM and the non-adiabatic pulsational code GraCo. We 
have applied the Frequency Ratio Method (FRM) and the Time Dependent Convec- 
tion theory (TDC) to estimate the mode identification, the Brunt- Vai'sala frequency 
integral and the mode instability, making the selection of the possible models. When 
the non-seismological constraints (i.e its position in the HR diagram) are used, the 
solar abundance models are discarded. This result contradicts one of the main hy- 
pothesis for explaining the A Bootis nature, namely the accretion/diffusion of gas by a 
star with solar abundance. Therefore, according to these results, a revision of this hy- 
pothesis is needed. The inclusion of accurate internal chemical mixing processes seems 
to be necessary to explain the peculiar abundances observed in the surface of stars 
with internal subsolar metallicities. The use of the asteroseismological constraints, like 
those provided by the FRM or the instability analysis, provides a very accurate de- 
termination of the physical characteristics of HR8799. However, a dependence of the 
results on the inclination angle i still remains. The determination of this angle, more 
accurate multicolour photometric observations, and high resolution spectroscopy can 
definitively fix the mass and metallicity of this star. 

Key words: stars: abundances - stars: chemically peculiar - stars: fundamental 
parameters (mass, metallicity) - stars: planetary systems - stars: variables: other - 
stars: individual: HR8799 



1 INTRODUCTION 

The A5 V star HR8799 (V342 Peg, HD 218396, HIP 114189) 
has been extensively studied in the last years. The 
first studies of this star were done in the context of 
asteroseismology. This technique is being developed as 
an efficient instru ment for the study of stellar interi- 
ors a nd evolution (|Aerts. Christensen-Dalsgaard. fc Kurt j 
I2OIOD . Through the comparison of the observed pulsa- 
tional modes with theoretical calculations, the general and 
internal character i stics of the star can be determined. 
ISchuster fc NissenI (|l986l ') firstly reported HR8799 as a pos- 
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sible SX Phoenicis type. IZerbi et al.l (| 19991 ) observed this 
star in a multisite multicolour photometric campaign, with 
Stromgren filters, and found three independent pulsational 
frequencies (/i = 1.9791 cd"^ (= 22.906 ^Hz), /2 = 
1.7268 cd"^ (= 19.986 AtHz), and /g = 1.6498 cd"^ (= 
19.095 A'Hz), units are cycles per day), making it one o f 
the 12 first 7 Doradus pulsators known l|Kave et al.lll999l ). 
This pulsating stellar group is composed of Main Sequence 

S MS) stars in the lo wer part of the classical instability strip 
Dupret et al. II2OO5'). Their pulsation modes have periods 
in the range [0. 5,3] days, that is, th ey are asymptotic g- 
mode pulsators. iGrav fc Kavd (|l999l ) obtained an optical 
spectrum of HR8799, and assigned an spectral type of kA5 
hFO mA5 V A Bootis, reporting an atmospheric metallic- 
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ity of [M/H]=-0.47. They also noted that HR8799 may be 
also a Vesa-typ e star, cha r acteri zed by a far IR excess due 
to a debris disk. 'Sadak and (|2006l ) developed a deep study of 
the metal abundances of this star, confirming its A Bootis 
nature (with surface chemical peculiariti es) . 

Two years later. [Marois et al] (|2008l ) reported the pres- 
ence of a planetary system around this star. It was the first 
detection of such a system carried out by direct imaging. 
This was the starting point of a set of approximately ten 
studies about this system during 2009, none of them from 
the asteroseismic point of view, even considering that the 
three detected pulsational frequencies can be used to better 
understand the host star, in particular its A Bootis nature 
and evolutionary s tatus. 

Discovered by [Morgan et akl j 19431 ). the A Bootis-type 
stars are non-magnetic, moderately-rotating, Population I 
stars with spectral types from late B - early A to F (dwarfs 
which show peculiarities in the morphology and abundance 
of the Fe-peak element lines. In particular, these lines are 
unusually weak considering their spectral types. Significant 
deficiencies in their abundances (up to 2d ex) are found, 
wher eas C, N, O and S have solar abundance ijPaunzen et al.l 
I2OO2I ). Different theories have attempted to explain the A 
Bootis nature from both observational (photometry, spec- 
troscopy) and theoretical investigations. It is not o ur aim to 
discu ss all of them here (for a interesting review see lPaunzenI 
I2OO3I ). Nevertheless, it may be worth describing the most 
accepted scenario relyi ng on the accretion of inter-stellar 
medium gas by the star (|Venn fc Lambertlll990l ). The accre- 
tion/diffusion scenario would explain the abundances found 
at the base of the outer convective zone of these stars, since 
convective layers are assumed to remain chemically homoge- 
neous. This accretion/diffusion scenario is based on the ac- 
cretion of inter-stellar medium by the star, and the mixing 
of these elements with those of the star due to diffusion and 
rotationally mixing processes. The accretion rate required 
to maint ain this situation is of the order of lO'^^-lO"^* M© 
per year ()Turcotte fc Charbonneciulll993l ). and, once the ac- 
cretion has ceased, the metal deficiencies should disappear 
in approximately 1 Myr due to diffusion and internal mix- 
ing processes. Thus, a possible interpretation is that A Boo- 
tis stars are young A-type stars (in a pre-main sequence or 
zero-age main sequence evolutionary stage) , still interacting 
with their primordia l clouds of gas and dust. Interestingly, 
iPaunzen et al] (|2002l ) found that most of the known A Boo- 
tis stars lie between the ZAMS (zero-age main sequence) 
and the TAMS (terminal-age main sequence, with ages of 
several hundreds million years). In this case the most likely 
scenario would be a MS star wi th solar abun dance pass- 
ing through an interstellar cloud (|Kamp fc Paunz cn 2002). 
However, the chemical mixing due to internal processes, such 
as rotationally-induced mixing, cannot be discarded as a 
possible explanation of the observed abundances. 

A Bootis stars and other types of objects in the same 
region of the HR diagram, such as 5 Scuti and 7 Doradus 
stars, are considered as particularly suitable for the aster- 
oseismological study of poorly known hydrodynamical pro- 
cesses occurring in stellar interiors, like the extent of the 
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convective core, mixing of chemical elements, redistribution 
of angular momentu m (jZahnl Il992l : iMathis fc ZahnI |2004| : 
iDecressin et al.ll2009l ). etc. A Bootis-type stars are also pul- 
sating stars. Therefore, asteroseismology can be used to ob- 
tain information about the internal structure of these ob- 
jects. Several works have been devoted t o study A Booti s 
stars with 5 Scuti pulsations, for instance iPaunzenI l|l998h : 
ICasas et alj (|2009l ). In addition, the combined use of space 
and ground-based obs ervations improves the potential of the 
modelling of the star (|Bruntt et al.ll2007l ). 

The present work aims at a comprehensive asteroseis- 
mic modelling of HR 8799, focusing on the discussion of the 
A Bootis nature of the star. More precisely, we want to an- 
swer whether the observed abundance is intrinsic to the 
star or an effect of surface processes. This answer drives 
the search for possible mechanisms explaining the observed 
abundances, that is, if there is accretion or not, together 
with internal chemical transport (rotationally-induced mix- 
ing, gravitational settling, radiative levitation, etc.). Up to 
now, only three A Bootis stars h ave been reported to b e 
7 Doradus pulsators: HD 218427 (|Rodn'guez et all l2006al ). 
HD 239276 (Rodriguez et al. 2006:3), and HR8799. As 7 Do- 
radus stars are asymptotic g-mode pulsators, they are very 
good candidates for testing the internal structure of the star, 
in particular its internal metallicity. Some of the most up- 
dated tools adapted for this purpose arc used: 1) the evolu- 
tionary code CESAM (.Morel fc Lebreton 2008). and 2) the 
pulsa tion code GraCo ( Mova et an[2004l : iMova and Garridd 
I2OO8I ). Using these tools we performed a massive numerical 
study of HR 8799 in an attempt of constraining physical and 
theoretical parameters. In this work we will follow the same 
scheme used for the study of RV Arietis, 29 Cygn is and 9 
Aurigae jCasas et al.ll2006l . l2009l : iMova et al.ll2006l ). 



2 OBSERVATIONS AND THE PROPERTIES 
OF HR8799 

The A Bootis star HR8799 was discovered to be a pulsating 
7 Doradus star bv lZerbi et al.l (| 19991 ). They observed it us- 
ing multicolour Stromgren photometry and found three fre- 
quencies with their corresponding amplitudes and phases in 
the different wavelenght bands. In Table 1 these frequencies, 
and the amplitudes and pha ses for the tw o high est ampli- 
tude frequencies, are shown. iGrav fc Kav"3 (|l999h obtained 
spectroscopic observations of this star, providing accurate 
values for the stellar luminosity, T^s, and log g (see Table 
2). They also established the A Boot i s and Vega-type nature 
of this star. In a study bv I Sadakand ()2006l ). a more detailed 
spectroscopic analysis of this star was done. An accurate de- 
termination of the individual atmospheric abundances was 
obtained, confi r ming the previous results. In a recent paper, 
ICuvpers et al.l l|2009l ) carried out a long-term photometric 
monitoring of a set of 7 Doradus stars, including HR8799. 
The effective te mperature reported in that work is similar 
to that given in lGrav fc Kavd (| 19991 ). 

An interesting additional property of HR8799 was the 
observation, by HST and Spitzer, of two well separated rings 
in its debris disk, locat ed at the inner and the outer limits of 
the p lanetary system |Chen et al.ll2009l : iReidemeister et al.l 
I2OO9I) . 
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Figure 1. Spectrum of HR 8799 in the interval 4000-5000 A 
(black) and, superimposed, a synthetic spectrum computed with 
Teff = 7430 K, logg = 4.35, [M/H]=-0.5 and Dsini=37.5 km s"! 
(red), based on model atmospheres by Kurucz. 



2.1 New spectroscopic observations 

We have obtained intermediate resolution, high signal-to- 
noise spectra in the range 3600-7800 A, with CAFOS (Calar 
Alto Faint Object Spectrograph) on the 2.2-m telescope at 
Calar Alto Observatory (CAHA, Almen'a, Spain). The ob- 
servations were taken on 9 November 2009. CAFOS was 
equipped with a CCD SITe detector of 2048x2048 pixels 
(pixel size 24 /xm) and the grisms Blue-100 and Green-100, 
giving a linear dispersion of ~ 88 A/mm (2 A/pixel). 

Spectroscopic observations in the range 3600-6800 A 
were also obtained on 12 November 2009 with the spectro- 
graph Albireo on the 1.5-m telescope at the Sierra Nevada 
Observatory (OSN, Granada, Spain). This instrument is 
equipped with a CCD2k detector of 2048x2048 pixels. The 
grating 1200 was used, providing a linear dispersion of ~ 67 
A/mm (1 A/pixel). 

In both cases the slit width was 1.5 arcsec. The usual 
bias, dark, dome flat-field and calibration lamp frames were 
taken. Standard procedures were used to process the data. 
The widths of the lines of the calibration lamps were used 
to obtain the instrumental responses of the spectrographs. 



2.2 Spectral fitting and properties 

The aim of the observations was to carry out a double check 
of some of the stellar parameters collected from the litera- 
ture. In Fig. 1 we show the CAFOS spectrum of HR 8799 
in the interval 4000-5000 A (containing the H/3, H7 and 
H5 Balmer lines) plotted as a black solid line, superim- 
posed by a synthetic spectrum computed using Tofi = 7430 
K, log5 = 4.35, [M/H]=-0.5 and i;sini=37.5 km s'^ The 
resolution of the synthetic spectrum has been degraded to 
match that of the observed one. The agreement between 
both spectra is quite remarkable, indicating that the stellar 
parameters used are reliable. 

The spectral synthesis has boon done using the AT- 
LAS9 and SYNTHE codes bv lKuru cz (1993). In this work we 
have used the GNU Linux v ersion of the codes available on- 
lin j^ llSbordone et al 1 12004|) . Model atmospheres provided 
bv ICastelU fc Kurucj (|2003l l are used as one of the inputs 
for the synthesis. 

Fig. 2 shows the position of HR8799 in the colour 
magnitude diagram (CMD). Evolutionary tracks are from 
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Claret! (Il995l). The m ain sequence has been taken from 
Philip fc Egre^ (Il980l) and the 5 Scuti instability strip 
from iRodrfgue z fc Brcgcr (2001) . Th e 7 D or instability 
strip is from Handler & Shob broo^ (|2002l ). Solid stars 

represent bona fide 7 Doradus objects observed by 

[Handler fc Shobbrookl l|2002l ), and solid triangles are bona 
fide 7 Doradus stars from the literature. The solid and open 
circles are prime and less probable 7 Doradus candidates re- 
spectively from those authors. We have derived photometri- 
cally My and {b — y)o using Str omgren data obtained from 
the Hauck-Mermilliod catalogue l|Hauck fc Mermilliod|[l998l . 
;solid star). The reddening correction was negligible. On 
the other hand, using the Hipparcos parallax (25.38 ± 0.85 
mas, i.e., a distance of 39.40 pc), an absolute magnitude 
My =5. 96 was calculated. The Hipparcos My is plotted as 
an open star. Other available parallaxes in the literature do 
not change significantly the absolute magnitude of this star 
(less than 0.1 L/Lq). 

The use of different scales of bolometric corrections 
yields significant larger luminosities (up to 0.3 Lq), com- 
promising the study presented here. We have used the 
most updated too l provided by Virtual Observatory, VOSA 
(jSavo et al.ll200^ to fit Kurucz atmospheric models to all 
the photometric observational data available in the litera- 
ture, avoiding the necessity for estimating the stellar lumi- 
nosity using bolometric corrections. Models with realistic 
metallicities, best fitting o bservations, p r ovide a luminos- 
ity similar to that given in iGrav fc Kav"3 l|l999l ) (about 0.1 
L/Lq larger). 

2.3 Rotation and the limits of our analysis 

An important observational uncertainty for the present 
study comes from the inclination angle i between the ro- 
tation axis of the star and the observer line of sight. Us- 
ing the rotation velocity v sini = 37.5 km s~^ given by 
iKave and Strassmeierl l|l998l ) , several wo rks in the literature 
tried to provide estimations for i, as in iReidemeister et al.l 
('2009'), where one of the detected frequencies have been 
used as signature of the rotation velocity, obtaining three 
estimations for the inclination angle, one per observed fre- 
quency. The amplitude and phase differences observed in 
different colour bands for these frequencies make this in- 
terpretation invalid. Other determinations of the inclination 
angle in the literature are related with the observation of 
the debris disk and/or the dynamics of the planetary sys- 
tem. iFabrvckv fc Murrav-Cla-^ (|2010l ) found a lower limit 
i > 20° based in the dynamical stability of th e system, as- 
sumin g the mass of the pl anets provided by iMarois et al.l 
(|2008l ). On the other hand. lSu et aD (|2009l ) studied the de- 
bris disk of HR 8799 and found that inclination angles larger 
than 25° must be ruled out. We have to take this latter re- 
sult with caution since we must take into account that the 
morphology inferred from the observations at 70 /xm is not 
a good constraint for the inclination angle of the star. That 
emission is dominated by a distribution of small grains and 
the geometry and the inclination angle we derive would cor- 
respond to those of disk-like distribution of such particles. 
That inclination angle may not coincide with the inclination 
angle of the star (|Bate. Lodato. fc Pringlj|2010l ). 

In the present study, we use the F reque ncy Ratio 
Method (FRM) developed in iMova et all (120051 ) (see Sec- 
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Table 1. Stromgren multicolour photometric observations. Phase differences and amplitude ratios of filters 
M, V and h related to filter y for the two frequencies with highest amplitude are also shown. 
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1.9791 {= 


22.906 AtHz) 


-0.5±3.6 


2.5±2.8 


1.8±2.8 


0.936±0.057 


1.470±0.071 


1.275±0.061 


1.7268 {= 


19.986 mHz) 


-7.0±7.5 


-0.8±5.7 


-1.4±5.7 


0.830±0.112 


1.316±0.137 


1.181±0.131 


1.6498 {= 
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Table 2. Physical chara cte ristics of HR 8799 taken 
iKave and Strassmeierl lll998l ') 3) Ivan LeeuwenI l|2007l 'l 



from: 1) iGrav fc Kav3 lll999l') 2) 



Teff (K) 


7430±75 


ref. 1 


log g (cms~'^) 


4.35±0.05 


ref. 1 


Mv 


2.98±0.08 


ref. 1 


R{Rq) 


1.34±0.05 


ref. 1 


L{Lq) 


4.92±0.41 


ref. 1 


V sin i (km s ^) 


37.5±2 


ref. 2 


TT (mas) 


25.38±0.85 


ref. 3 





(b-y). 
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Figure 3. Rotation and angular velocities of HR8799 as a func- 
tion of the inclination angle i. The red region (i < 18°) represents 
the zone where the star can hardly be a A Bootis star. The green 
region (18° < i < 36°) represents the zone where our proce- 
dure cannot be applied accurately. The horizontal lines are the 
observed frequencies of HR8799. 



Figure 2. Colour magnitude diagram with absolute visual mag- 
nitudes My plotted agai nst dereddened Stromgren (b —y)o colour. 
Small symbols are from [Handler &: Shobbrookl ||2002| '): solid stars 
represent bona fide 7 Doradus stars observed in their work, and 
solid triangles are bona fide 7 Doradus stars from the literature. 
The small solid and open circles are prime and less probable 
7 Doradus candidates respectively. The position of HR8799 is 
determined by using its Stro mgren magnitud es taken from the 
Hauck-Mermilliod catalogue ijHauck &: Mermi lliod 1998) and its 
absolute visual magnitude. Two values were derived for the lat- 
ter, one der ived from the St rom gren photometry using the cali- 
brations bv lCrawfordI l ll975l) andlCrawfor 3 1119791) . implemented 
in the program TempLogG llRoeerslll99'5ir and another one using 
the Hipparcos parallax. These two determinations are given in 
the plot by the larger solid and open star symbols respectively. 



tion 3.2) for the asteroseismological modelling of this star. 
This procedure has a limited rang e of validity in stell ar ro- 
tation velocity, as demonstrated in lSuarez et al.l(|2005l ). The 
rotation velocities used in this work and the pulsational fre- 
quencies are around 0.1 and 0.25 of the break-up rotation ve- 
locity respectively. Therefore, the use of a second order per- 
turbation theory for the calculat ion of the frequencies with 
rotation is accurate enough (Ba llot et al. Ilioog). That sec- 
ond order theory is used for the study done bv lSuarez et 
(j2005j). Nevertheless, as the observed frequencies are larger 
but close to twice the angular rotation velocities used in this 
work, the use of the traditional approximation to study the 
accuracy of the asymptotic exp ression in the pre s ence of ro- 
tation would be rec ommended (|Lee fc Saidll987l :l Townsendl 
I2OOOI : iMathidliooi ). This study needs to be carried out to 
determ ine more prec i sely t he reliability of the FRM. 

Ini uarez et all (|2005l ). the authors found a maximum 
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rotation velocity of around 60kms~^ for the correct ap- 
plication of th e FRM, for a standard 7 Dorad us star. On 
the other hand, iTurcotte fc Charbonneaul (!l993') found that 
meridional circulation cannot destroy the accretion pattern 
(in the accretion/diffusion scenario) for a rotation velocity 
below 125kms~^. These values provide two lower limits for 
the inclination angle, one limited by the A Bootis nature and 
another one to ensure the applicability of the FRM. 

In Figure 3, the rotation velocity as a function of i for 
a stellar radius of ii = 1.30 7?0 is shown. For the different 
stellar radii determined for this star the results are similar. 
In this figure we see in red the region with i < 18°, limit 
imposed by the A Bootis nature of HR8799. In green it is 
the region with 18° < i < 36° , the Umit for the FRM to be 
applicable. We want to note that if the inclination angle lies 
in the green zone, our analysis would be inaccurate. Finally, 
the observed pulsational frequencies are displayed with hor- 
izontal lines. In this figure we see how those frequencies are, 
at least, twice the rotation frequencies used in this study. 



3 TOOLS AND METHODS 

In this section, we give a brief description of the tools used 
for the present study. 



3.1 Stellar equilibrium models 

The stellar equilibrium models were computed using the evo- 
lutionary code CESAM, with a mesh grid (B-splines basis) of 
2000 points. First-order effects of rotation on the equilibrium 
models were considered by subtracting the spherically aver- 
aged contribution of the centrifugal acceleration to the grav- 
ity of the model, QcB ~ g — Ac{r), where g corresponds to the 
local gravity, and Ac{r) represents the radial component of 
the centrifugal acceleration. This spherically averaged com- 
ponent of the centrifugal acceleration does not change the 
order of the hydrostatic equilibrium equa tions. Such mod- 
els are called 'pseudo-rotating' models (see lSoufi et al. 1998: 
ISuarez et al.l [2007.1 . During evolution, models are assumed 
to rotate as a rigid body, and their total angular momen- 
tum is conserved. Although the non-spherical components of 
the centrifugal acceleration were not considered, they were 
included as a perturbation in the computation of the oscil- 
lations. 

Standard physical inputs for 7 Doradus 
stars are used, i.e., the CEFF eq uation of state 
ijChristensen-Dalsgaard fc Dapper] Il992l ) . The opacity ta- 
bles w ere taken from the OPAL package (|lglesias fc Rogeri 
Il996l ). complemented a t low temperatures (T ^ 10"^ K) by 
the tables provided bv [Alexander fc FergusonI (|l994l '). The 
atmosphere was calculated following two approaches: a grey 
atmosphere (Eddington T{t) law) when the equilibrium 
models were used to compute adiabatic oscillations, and 
the Kurucz model atmospheres for the computation of 
non-adiabati c quantities. The abund ance mixture used is 
that given in lGrevesse fc Noelj (119931 ). 

The main approximation taken in the models with a 
possible infiuence in the results of this study is the lack of 
updated internal chemical transport mechanisms in the equi- 
librium models, such as rotation-induced mixing or gravita- 
tional waves. The inclusion of these mechanisms is the next 



step to be taken, considering the conclusions obtained in the 
present work. 



3.2 Oscillation computations 

The seismic models were completed by computing for each 
equilibrium model its corresponding oscillation spectrum 

and non-adi abatic observable s. 

GraCo (|Mova et al.ll2004lMov a and Garridd[2008h pro- 
vides non-adiabatic quantities related to pulsation and in- 
cludes the convection - pulsation interaction using th e Time 
Dependent Convection theory l|Dupret et al.l |2005| ). This 
theory improves the "frozen convection" approximation usu- 
ally implemented in most of the codes. With this we can 
study the instability of the modes. It also provides accurate 
variations of the effective temperature, effective gravity and 
phase-lags of the modes. These quantities are the base of 
the mode identification using multicolour photometry (see 
Section 3.5). The adiabatic solutions of this code have been 
used as reference for the ESTA works (Evolution and Aster- 
oseismic Tools Activities, iLebreton et al. I l2008l : iMova" et al.l 
l2008h . 



3.3 The Frequency Ratio Method 

The Frequency Ratio Method or FRM l|Mova et al.l |2005| ) 
takes advantage of the analytical expressi on of the fre- 
quencies given by the asymptotic theory (|Tassoull Il980l : 
b movers fc MoTO H2007h . 7 Doradus stars and other slow pul- 
sators present frequencies in this g-mode asymptotic region. 
This makes it possible to obtain information of the radial or- 
der n of the modes and the Brunt- Vaisala frequency integral 
of the model, through the ratios of the observed frequencies. 
With at least three frequencies the method can be carried 
out. 

The application of the method to the three observed fre- 
quencies provides chains of {ni, £i,n2, i2,n3, (.3, lobs ) , where 
Hi is the radial order of mode i, £i is its spherical order, 
and lobs is the new observable introduced by the FRM, the 
observed Brunt- Vaisala frequency integral, the same for the 
three frequencies since it is a global property of the stellar 
model. 

iMova eTaLl l|2005h showed that the real solution is part 
of the set of possible chains provided by the FRM. In ad- 
dition, in this work we have relaxed the requirement of all 
the modes having the same £. We allow £ to take values 1 or 
2, the most probable ones for modes detected with ground- 
bas ed photometry for 7 Doradus stars. 

ISuarez et al. (|2005l ) showed the application limits of the 
FRM due to stellar rotation. With the assumption of all the 
frequencies having the same azimuthal order m, the main 
conclusions of that paper were: 1) the rotation limit for a 
standard 7 Doradus star for the perturbative approxima- 
tion to be valid is around 60 kms~^, 2) if that assumption 
is fulfilled, the real solution is still part of the list of possible 
solutions provided by the FRM, otherwise the FRM does 
not give any solution, and 3) the accuracy of the asymptotic 
expressions decreases with the rotation velocity, increasing 
the number of possible solutions. Nevertheless, these con- 
clusions and application limits of the FRM must be revised 
under the traditional approximation (see section 2.3). 
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3.4 Instability analysis 

The energy balance of the modes during a complete period 
is one of the most useful information we can use. Analysing 
the energy balance of a pulsational mode with its surround- 
ings in a pulsation cycle we can infer if the mode is sta- 
ble (losses energy) or over-stable (gains energy). To do this 
study accurately, the main physical processes having influ- 
ence on this energy bala nce must be implemente d in the 
codes. iGuzik et al'l (| 200dh andl Pupret et al.1 l|200d ). showed 
that the position of the base of the the convective envelope is 
the key point driving 7 Doradus modes, since they are origi- 
nated by the blocking of the radiative flux at this point. This 
can explain the location of the blue edge of the 7 Doradus 
instabili ty strip. For these s t ars, the Time Dep endent Con- 
vection (jPupret et al.ll2005l ; iGrigahcene et al] 2005. TDC) 
is the only set of equations predicting unstable modes in the 
observed frequency ranges. This development has provided 
the first theoretical predictions of instability strips for these 
stars. 
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Figure 4. Position on the HR diagram of the models fulfilling 
the spectroscopic observations. The label "nosis" at the top of 
the figure means that no seismological constraints are used. 



3.5 Multicolour photometric analysis 

In the following, multicolour photometry is used to pro- 
vide constraints on the physical characteristics of the star 
through additional information on the degree I of the spher- 
ical harmonic associated to each observed pulsational fre- 
quency. 

The linear approximation to non-radi al flux varia- 
tions of a pulsating star was first derived by Dziembowskil 
lll977h, and later reform ulated by iBalona fc Stobid (|l979f ) 
and lWatsonI (|l988l ). Then. lGarrido et al.l l|l990l ') showed that 
the V and y Stromgren bands can be used for discriminat- 
ing the degree i. The comparison of the numerical solutions 
with the observa tions is based on n on-adiabatic calculations 
(more details in iMova et al.l |2004| ) . In particular, pulsation 
is non-adiabatic in the outer convective zone of these stars, 
where the thermal relaxation time is of the order of the pul- 
sation period. Accurate determination of the eigenfunctions 
in these layers, therefore, requires the use of a non-adiabatic 
description which includes the already mentioned TDC. This 
procedure makes it possible to relate multicolour photomet- 
ric observables with such eigenfunctions, allowing therefore 
a direct constraining on some unknown physical parameters 
through the direct comparison with observations. 

All these tools can be combined to provide an accurate 
determination of the general physical properties of IIR8799, 
as described in the next section. 



4 ANALYSIS OF THE STAR 

The analysis of this star has been done in different steps, 
imposing the observational constraints sequentially. The 
sequence is: 1) Fulfill the physical characteristics given 
by spectroscopy (Grav & Kave 1999), 2) fulfill the Brunt- 
Vaisala frequency integral constraints provided by the FRM, 

3) predict the observed frequencies to be over-stable, and 

4) identify the spherical degree I of the mode for the fre- 
quency with the highest amplitude by using the observed 
multicolour photometry. 



4.1 Position in the HR diagram 

The first step of our study was to search for models in a 
dense grid fulfilling the observed physical characteristics of 
the star. As it has been explained in Section 2, we have used 
the effective temperature, luminosity, and gravity (with their 
errors) given in Table 2. The internal metallicity has been 
regarded as a free parameter to test the A Bootis nature 
of the star and to obtain the internal metallicity using as- 
teroscismology. The mass, estimated in 1.47 ± 0.30 M© by 
iGrav fc Kavd (| 199^1 ). has been also regarded as a free pa- 
rameter since it has not been directly determined. 

We have developed a grid of CESAM equilibrium mod- 
els varying the mass (in the range [1.25, 2.10] M© with steps 
of O.OIM0), the metallicity (with values [M/H]=0.08, -0.12, 
—0.32, and —0.52), the Mixing- Length parameter qmlt (val- 
ues 0.5, 1, and 1.5), and the overshooting (values 0.1, 0.2, 
and 0.3). 

Figure 4 shows in an HR diagram the models in our 
grid fulfilling the observations. In Table 3, the range of free 
parameters of acceptable models is shown. With this first set 
of observational constraints we have reduced the number 
of acceptable models from around 600, 000 to 1975. This 
number is still large, but the first main result of this study 
is that no models with metallicity larger than solar are in the 
set of acceptable models. To study the reason for this, a new 
set of models with solar metallicity has been computed. Most 
of these tracks present models fulfilling all the constraints 
except the luminosity. In Fig. 5 we present the luminosity as 
a function of the effective temperature for tracks of different 
metallicities and masses. All the solar metallicity models 
present luminosities larger than the observed one for the 
effective temperatures adopted (represented as a box in the 
figure) . 

We conclude, then, that an internal solar metallicity 
must be discarded for this star. A luminosity 0.2 L/ Lq larger 
or a lower effective temperature can make solar metallicity 
models marginally acceptable. These would be young models 
located at the ZAMS. 

The next goal of this study is to obtain an internal 
metallicity as accurate as possible with the help of astero- 
seismology. 
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Table 3. Acceptable models depending on the physical constraints. Each row shows the results obtained 
when a new constraint is added with respect to the previous one. The bottom row shows the results obtained 
when all the constraints are imposed to the less favorable case from the FRM point of view (see text for 
details). The mass and metallicity determinations are linked, that is, the first mass range is linked with the 
first metallicity shown, etc. 



Constraint Mass Metallicity c*mlt Overshooting 

(in M0) (in [M/H]) 



HR position 


[1.25, 1.27], [1.32, 1.35], [1.40,1. 48] 


-0.52, 


-0.32, 


-0.12 


0.5, 1.0, 1.5 


0.1, 0.2, 0.3 


Previous + FRM 


[1.25,1. 27], [1.32, 1.34], [1.40,1. 48] 


-0.52, 


-0.32, 


-0.12 


0.5, 1.0, 1.5 


0.1, 0.2 


Previous + Inst 


[1.32, 1.33], [1.44,1. 45] 




-0.32, 


-0.12 


1.5 


0.1, 0.2 


Previous + colours 


1.32 






-0.32 


1.5 


0.1, 0.2 


Total in less favorable case 


[1.32, 1.33], [1.44,1.45] 




-0.32, 


-0.12 


1.5 


0.1, 0.2, 0.3 
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M=1.7, [Fe/H]=0 
M=1.60, [Fe/H]=0 
M=1.50, [Fe/H]=0 
M=1.48 , [Fe/H]=-0.12 
M=1.43 , |Fe/H]=-0.12 
M=1.45 , [Fe/H]=-0.32 
M=1.40,|Fe/H|=-0.32 
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Figure 5. Effective temperature vs. stellar luminosity of some 
evolutionary tracks (from PMS to ZAMS) with different masses 
and metallicities, as described in the labels. The box represent 
the location of valid parameters for HR8799 



4.2 Frequency Ratio Method 

The observation of three pulsational frequencies in the 
asymptotic g-mode region makes the FRM suitable for the 
seismological study of this star. First we must verify if the 
observed frequencies can be part of the same rotation split- 
ting (making the FRM not valid). This verification is made 
by taking into account that the largest difference among the 
observed frequencies is 3.8 /iHz. On the other hand, the first 
order approximation to the rotation splitting provides a dif- 
ference of the order of the rotation velocity, weighted by a 
constant usually l arger than 0.5 for £ = 1 g-modes and 0.8 
for £ = 2 g-modes (|Chlebowskl 19781 ). In Figure 3, the small- 
est rotation velocities, obtained for i « 90° , are larger than 
6.52 /iHz. Therefore, the possibility of the two more sepa- 
rate observed frequencies being part of the same splitting is 
unlikely, since we would need the combination of a Ledoux 
constant around 0.6 (only possible for some i — 1 asymp- 
totic g-modes) and, an inclination angle larger than 70°, a 
very unlikely value for this star taking into account the im- 
age obtained of this system, suggesting that the inclination 
angles are not close to equator-on rotation. This means that 
the first applicability condition of the FRM is fulfilled. 

The second verification of the appl icability of th e FRM 
is done on the rotation velocity itself. ISuarez et al.l (|2005l ) 
showed that, as the analytical asymptotic expression is de- 
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Figure 6. Effective temperature vs. Observed Brunt- Vaisala in- 
tegral (lobs) obtained with the models fulfilling the spectroscopic 
constraints (upper panel, no seismological constraints used "no 
sis") and those restricted by the FRM (bottom panel, seismolog- 
ical constraints "sis"). 



veloped assuming no rotation, the uncertainty in the appli- 
cation of the expressions grows with the rotation velocity, 
until a limit where the method cannot be used (as explained 
in section 3.3). 

One of the main assumptions of the FRM is the need 
that all the frequencies used must have the same azimuthal 
order m. Currently there is no method that can provide an 
univocal determination for the value of m, the only one pro- 
viding some information being optical high-resolution spec- 
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troscopy. Therefore, we cannot verify this assumption obser- 
vationally, although there are some ind ications of th e fulfil l- 
ing of this requirement: 1) The work of ISuarez et al.l l|2005h . 
where the case of modes with different m is studied. Their 
conclusion is that the procedure is self-consistent. The FRM 
gives, in that case, no-solutions, but not a wrong solution. 
This conclusion can be, nevertheless, conditioned by the use 
of the perturbative approximation. 2) The visibility of the 
modes with different m in a rotating star is described in 
iGizon fc Solankil (|2003l '). For the inclination angles studied 
in this work, the visibility of the modes favours those with 
m = ±1 in the case of the only mode identified for this star 
(/i having 1 = 2, see next sections). Therefore, the only pos- 
sibility for the observed frequencies to have different m is to 
have m=l and m = —1. Taking into account the previous 
study of the value of the rotational splitting, this means that 
the respective m=0 are well separated and models can only 
marginally over- stabilize both modes at the same time. In 
the unlikely case of quasi equator-on inclination angles, the 
visibility of the modes favours those with m = ±2, resulting 
in a larger separation of the corresponding m = modes 
and an even more unlikely over-stabilization of both modes 
at the same time. The only possible scenario invalidating 
this assumption is, then, a /2 having 1 = 1 and m=0. 

In this work we present two studies for two possible 
rotation velocities. A first one for a moderate rotation ve- 
locity from the point of view of the FRM, around 45 kms~^, 
(this value has been chosen just for illustration purposes be- 
cause it clearly shows how the p rocedure used here works). 
iFabrvckv fc Murrav-Clavl (|2010l i found that this rotation ve- 
locity cannot be excluded for HR8799, acording to the dy- 
namical study of the system. In that case the complete pro- 
cedure gives very restrictive results regarding the number 
of accepted models. At the end of the paper (Section 4.5) 
we show the results obtained for the largest rotation veloc- 
ity admitted by the FRM (around 60 kms~^), where this 
procedure is less restrictive. 

In the case of a rotation velocity of 45 kms~^ we have 
obtained a set of 10 chains, each one with a fixed set of 
mode identifications and Brunt- Vaisala frequency integral. 
Fig. 6 shows the set of models fulfilling the additional con- 
straint of lobs (bottom panel) as compared with those only 
fulfilling the spectroscopic observations (upper panel). The 
number of rejected models, for this moderate rotation ve- 
locity, is around 70%, but a large number of models is still 
accepted (see Table 3), and only those with overshooting 0.3 
are rejected. This dependence between the Brunt- Vaisala 
frequency and the overshooting parameter is well known 
l|Miglio et al.|[2008h and is part of the valuable information 
provided by the FRM. 



4.3 Instability analysis. 

The next step is to verify if the observed modes are predicted 
to be over-stable for the remaining models, as explained in 
Section 3.4. Unfortunately, the TDC is not a complete de- 
scription of the pulsation - convection interaction, and the 
energy balance provided by this study is not totally accu- 
rate. Therefore, we have adopted a conservative criterion 
for the acceptance of a model: if the modes in the observed 
frequency range present a growth rate positive or negative 
but very close to zero, they are accepted. With this crite- 
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Figure 7. HR, diagram position of the models fulfilling the spec- 
troscopic observations -|- FRM -|- instability analysis. 



rion, some incorrect models will be considered as acceptable, 
but we prefer this inconvenience instead of rejecting correct 
models. 

This study is robust in selecting models, keeping only 
10% of those filtered in the previous step (see Fig. 7 com- 
pared with Fig. 4). As the energy balance is a phenomenon 
mainly thermodynamic, the accepted models are located 
in a certain Tofj range, in the coolest part of the HR dia- 
gram. This instability study has also rejected all the models 
with metallicity —0.52, and Qmlt 0.5 and 1.0 (as expected). 
The remaining valid masses are in the ranges [1.32,1.33] and 
[1.44,1.45] M0. 



4.4 Mode identification with multicolour 
photometry. 

The last step of our study is to use the mode identifica- 
tion from the mult icolour photometry obtained for this star 
(Zer bi et al. 1 1 19991 ). As explained in Section 3.5, the non- 
adiabatic code, with the TDC implemented, makes it pos- 
sible to obtain accurate non-adiabatic observables for these 
stars. With these observables we can carry out the identifi- 
cation of the spherical degree I of the observed modes com- 
paring the theoretical predictions for different modes with 
observations. In Fig. 8, this comparison, in the phase differ- 
ences vs. amplitude ratio diagram, is done for a representa- 
tive model. As all the acceptable models at this stage are in 
a narrow range of temperatures and masses, the rest of the 
models provide similar results. 

In Fig. 8, the observational errors are displayed as 
boxes, and only the comparison between bands h and y is 
made. The comparison of bands v and u with y produce the 
same type of behaviour. It can be easily checked that the 
mode /2 cannot be identified due to its large observational 
errors, but /i is always identified as = 2. This is a strong 
constraint to the models, since the FRM predicts fixed mode 
identifications for each particular lobs ■ Therefore, we have to 
select only those models for which /i is predicted to have 
1 = 2. These models are shown in Fig. 9. 

The number of final models fulfilling all the observa- 
tional constraints is drastically decreased, remaining only 
22 models of the initial 1975 models selected not using 
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Mode identification using multicolour photometry 
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Figure 8. Phase differences vs. amplitude ratio between bands b 
and y of the two observed frequencies with highest amplitude 
(boxes) as compared with theoretical predictions for different 
spherical degrees I (points). 
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Figure 10. IfR, diagram position of the models fulfilling the spec- 
troscopic observations + FRM + instability analysis + multi- 
colour photometric mode identification, for the less favourable 
case in the FRM. 
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Figure 9. IfR diagram position of the models fulfilling the spec- 
troscopic observations -|- FRM -|- instability analysis + multi- 
colour photometric mode identification. 



asteroseismology. These 22 models have [M/H] =— 0.32 ± 
0.1 (higher than that given by iGrav fc Kav3 (j 19991 ) 
[M/H] =— 0.47 ± 0.1, but consistent with the error bars), 
M=1.32 ± O.OIM0, aMLT = 1.5 ± 0.3, and overshooting 
0.15 ±0.10. 



4.5 The less favourable case 

We have already mentioned that the rotation velocity is one 
of the main uncertainties for this procedure. All the previous 
studies have been developed assuming a rotation velocity in 
the middle of the applicability range of the method, that is, 
with a moderate rotation velocity, or an inclination angle 
i — 50° -60°. This may not be the case. Therefore we have 
also repeated the same procedure for the less favorable case: 
a rotation velocity around 60 kms^^, which implies an incli- 
nation angle around i « 36°. In this case, the most affected 
part of the procedure is the FRM, increasing the errors and, 
as a consequence, the number of acceptable models and lobs. 
The final result is not as detailed as for the moderate rotat- 



ing case, and the final set of models fulfilling all the con- 
straints is displayed in Fig. 10. The number of acceptable 
models is 66, with [M/H]=-0.32 and -0.12, M=[1.32,1.33] 
and [1.44,1.45] M0, omlt = 1.5, and overshooting= 0.1, 0.2, 
0.3. Analysing the differences with the v^ot ~ 45kms~^ case 
we have noted that the main one is the inclusion of models 
with overshooting 0.3, making it possible to fulfill the con- 
straints the models with [M/H]=-0.12, and M= [1.44, 1.45] 

Mq. 



5 CONCLUSIONS AND PERSPECTIVES 

In this work, the first comprehensive asteroseismologic study 
of the planetary system host HR8799, a A Bootis star pre- 
senting 7 Doradus pulsations has been carried out. This as- 
teroseismic work is specially important for the determina- 
tion of the internal abundances of this kind of stars, a previ- 
ous step to understand the physical mechanism responsible 
for the surface chemical peculiarities of the A Bootis group. 
The asymptotic g modes observed in HR8799 make it pos- 
sible to deeply test its internal metallicity, something hard 
to do with 5 Scuti pulsators, for example. 

For the asteroseismologic study, we have used part of 
the most updated codes and physics for describing 7 Do- 
radus stars: CESAM and GraCo codes, TDC theory, FRM, 
and mode identification using Stromgren multicolour pho- 
tometry. We have constructed a dense grid of equilibrium 
models, and carried out a first selection using the available 
data (TcfT, log g, and luminosity), leaving the metallicity and 
mass as free parameters in order to test the internal metal- 
licity of the star. 

The first selection shows that there are no models with 
solar metallicity fulfilling the observations, since the stel- 
lar luminosity derived from the observations is smaller than 
any of the possible luminosities of models with solar metal- 
licity. This contradicts the main assumption of the theories 
explaining the A Bootis nature, i.e. that these stars have so- 
lar metallicity, whereas the observed abundances are due to 
surface phenomena. 

The asteroseismic study to fix the metallicity of the star 
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is done using 1) the FRM, to estimate possible mode iden- 
tifications and Brunt- Vaisala integrals, 2) energy balance to 
determine the stability of the modes, selecting models in a 
certain temperature range and qmlt and, 3) mode identi- 
fication with multicolour photometric observables, selecting 
only models for which FRM predicts the mode with the 
highest amplitude to have £ = 2. This is the first time that 
such a study has been done for a plan et-hosting star, bu t 
the general procedure was proposed in iMova et al.l (|2008l ) . 
All these steps provide a small group of acceptable models 
depending on the rotation velocity. For a moderate rotation 
velocity (implying i around 50°) the mass and the metal- 
licity are very precise determined ([M/H]=— 0.32 ± O.f and 
M=1.32 ± O.OfM©). For a rotation velocity at the limit of 
the validity of the FRM (around i ~ 36°), the method is less 
selective and provides some possible values of the metallicity 
and mass ([M/H]=-0.32 and -0.f2±0.f , and M=[f .32,1.33] 
and [1.44,1.45] Mq respectively). We point out that inclina- 
tion angles lower than 18° cannot be possible to preserve the 
A Bootis nature of the star. Nevertheless, this limit must be 
also verified with more accurate observations. 

These results have some consequences. The first one 
is related to the most accepted theory of the A Bootis 
nature, namely the acc retion /diffusion scenario given by 
IVenn fc LambertI ([l990). The result of this work, discard- 
ing solar metallicity as the internal metallicity of the star, 
together with the fact that some A Bootis stars hav e debris 
disks not connected with the star (|Chen et al.l 20091) . make s 
this explanation unlikely, but not neglig ibl e JSu et al]|2009l) 
In an y case, the equation developed by IVenn fc Lambert 
|l99(7) to describe the individual abundances as a sum of 
internal abundance plus that of accreted material, must be 
corrected to include a possible non-solar internal abundance 
as following 



e(m) = (1 - /)e,(m) + /eisM(m) 



(1) 



where / is a mixture factor, and e, e« and eisM are the abun- 
dances observed, internal and coming from the interstellar 
medium of each chemical element, respectively. Therefore, 
the study of internal chemical mixing processes seems to be 
the key to explain the A Bootis nature, at least for HR8799, 
as the solar abundances for C, N, O and S observed on its 
surface have still to be explained. 

On the other hand, we have seen that the asteroseismo- 
logical study needs an accurate determination of the rotation 
velocity, not just the projected rotation velocity, of the star, 
for a very accurate determination of the stellar characteris- 
tics. Therefore, a consequence of this study is the need for a 
precise determination of the inclination angle i, of the mul- 
ticolour photometric amplitudes and phases of /2, and some 
information of m values through time-series if high resolu- 
tion spectroscopy. These determinations would help to carry 
out a definitive selection of the models. 
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